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Q-i, It has been shown that hybrid inflation may end with the formation of non- 

topological solitons of inflaton field. As a first step towards a fully realistic 
picture of the post-inflation era and reheating in super symmetric hybrid inflation 

X 

models, we study the classical scalar field equations of a supersymmetric hybrid 
inflation model using a semi-analytical ansatz for the spatial dependence of the 
fields. Using the minimal D-term inflation model as an example, the inflaton 
field is evolved using the full 1-loop effective potential from the slow-rolling era 
to the U(1)fi symmetry-breaking phase transition. Spatial perturbations of the 
inflaton corresponding to quantum fluctuations are introduced for the case where 
there is spatially coherent U(1)fi symmetry breaking. The maximal growth of 
the dominant perturbation is found to depend only on the ratio of superpotential 
coupling A to the gauge coupling g. The inflaton condensate fragments to non- 
topological solitons for X/g > 0.09. Possible consequences of non-topological 
soliton formation in fully realistic SUSY hybrid inflation models are discussed. 
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1 Introduction 



An important event in the history of the Universe is the process by which it became 
hot after inflation, known as reheating. A full understanding of this epoch is essential 
for a full understanding of cosmology. It is also important to have a clear picture of 
the post-inflation cosmological enviroment against which ideas in particle physics and 
cosmology can be judged. The conventional view of reheating following inflation has 
been based on the perturbative decay of a spatially homogeneous Bose condensate of 
inflaton scalars [I]. This picture, however, has been challanged in recent years. It 
was firstly observed that the inflaton condensate might decay non-perturbatively via 
parametric resonance [21 El H] ■ More recently, for some inflation models, in particu- 
lar hybrid inflation jS], it was shown that the inflaton sector Bose condensate can be 
unstable with respect to spatial perturbations of the scalar fields jHlIZj. The growth 
of the spatial perturbations could either be slow, in which case a coherently oscillat- 
ing condensate would form prior to spatial perturbations going non-linear ('inflaton 
condensate fragmentation') [6J, or rapid, in which case the spatial perturbations dom- 
inate the total energy density before the inflaton can enter into coherent oscillations 
('tachyonic preheating') [7\ In the case of inflaton condensate fragmentation it was 
observed that the final state of the process is likely to be, at least for some range of pa- 
rameters of the model, the formation of non-topological soliton-like objects ('inflaton 
condensate lumps') jUllE]- The formation of such objects in hybrid inflation models 
was subsequently confirmed in lattice simulations [UJ. Recently it has been shown 
that inflaton condensate fragmentation can also occur in single field chaotic inflation 
models [TP] . 

However, there has been to date no detailed study of condensate fragmentation 
and non-topological soliton formation in a fully realistic hybrid inflation model. Of 
particular interest is the case of supersymmetric (SUSY) hybrid inflation models, which 
are prime candidates for natural inflation models in the context of SUSY ^T] ^3 • (The 
alternative, SUSY chaotic inflation, is disfavoured by the initial inflaton expectation 
value required, larger than the natural scale of supergravity (SUGRA) corrections, 
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M = Mpi/\^8n [13J.) In the following we will focus on the case of the minimal D- 
term hybrid inflation model as a specific example. In D-term inflation models the 
energy density is provided by a Fayet-Illiopoulos D-term associated with a U(1)fi 
gauge symmetry ^2] ■ D-term inflation is particularly favoured as it provides a natural 
solution to the //-problem of SUGRA inflation models [T^j . 

Hybrid inflation models in general are based on two scalar fields, an inflaton, S, 
and a second field, $, which undergoes a symmetry-breaking phase transition which 
terminates inflation. (We will refer to these as the 'inflaton sector' fields.) As well 
as the homogeneous scalar fields considered in conventional D-term inflation models, 
there will also be quantum fluctuations leading to spatial perturbations of the scalar 
fields jl]. In order to establish the initial conditions for the growth of spatial perturba- 
tions following inflation, it is essential to follow the evolution of the inflaton from the 
slow-rolling epoch during inflation until the phase transition which terminates hybrid 
inflation. The evolution of the inflaton prior to the phase transition in SUSY hybrid 
inflation models is dominated by the 1-loop effective potential. Therefore in the fol- 
lowing we will consider the evolution of the inflaton sector fields in D-term inflation 
including the full 1-loop radiative correction. 

Spatial perturbations of the inflaton sector fields arise initially as a result of quan- 
tum fluctuations. A method to study the evolution of these fluctuations into the 
classical regime has been developed in [T5] . So long as the final state of the evolution 
of the scalar fields is semi-classical, the evolution may be studied by solving the clas- 
sical equations of motion but with a classical probability distribution for the initial 
conditions, determined by the initial quantum state of the field. (Following jT^j, we 
will refer to this as the 'equivalent classical stochastic field' (ECSF) method in the 
following.) 

In this paper we study the evolution of spatial perturbations in D-term inflation 
using a semi-analytical model with simplified initial conditions for the spatial pertur- 
bations. This will allow us to gain a physical understanding of the process of spatial 
perturbation growth whilst providing an example against which future numerical sim- 
ulations may be tested. We will argue that the growth of perturbations of the inflaton 
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in the semi-analytical model represents the maximum possible growth in realistic mod- 
els. In addition, we will show that that the formation of cosmic strings is likely to be 
rare on scales smaller than the horizon and so is unlikely to play a significant role in 
the subsequent evolution of the energy density, in which case the model we consider 
should give a good estimate of the growth of perturbations in the realistic case. 

As in all previous studies, we will consider the evolution of a real inflaton field. 
However, in fully realistic models of SUSY inflation, the inflaton is a complex field. As 
such, none of the existing analyses of spatial perturbation growth in hybrid inflation 
models fully describes the SUSY case. One result of the complex nature of the inflaton 
field is that the non-topological solitons which form at the end of hybrid inflation 
are likely to be Q-balls PHI |2Hj- We will discuss how the analysis presented here 
relates to Q-ball formation and discuss some possible cosmological consequences of 
non-topological soliton formation in SUSY hybrid inflation models. 

The paper is organised as follows. In Section 2 we discuss the D-term inflation 
model and the full 1-loop effective potential. In Section 3 we discuss the initial con- 
ditions used for the spatial perturbations. In Section 4 we discuss the semi-analytical 
ansatz used to model the growth of spatial perturbations. In Section 5 we discuss the 
evolution of perturbations of the inflaton sector fields in the semi-analytical model. In 
Section 6 we discuss some possible consequences of non-topological soliton formation 
in realistic SUSY hybrid inflation models. In Section 7 we present our conclusions. 
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2 Equations of Motion and 1-loop Effective Poten- 
tial 



In the following we will consider the minimal SUSY D-term hybrid inflation model 1 
[T2] . The superpotential of the minimal D-term inflation model is 

W = AS$+$_ . (1) 

The tree-level scalar potential is then 

V = A 2 |S| 2 (|$+| 2 + |$_| 2 ) +A 2 |$+| 2 |$_| 2 + |-(|<I> + | 2 - |$_| 2 + £) 2 , (2) 

where S is the inflaton, Q± are fields with charges ±1 with respect to a Fayet-Illiopoulos 
U(1)fi gauge symmetry, £ > is the Fayet-Illiopoulos gauge term (£ x / 2 « 8.5 x 
10 15 GeV ^21) and g is the U(1)fi gauge coupling. Since the minimum of the potential 
as a function of S is generally at \$+\ = and since |$ + | = at the end of inflation, 
we will assume this value throughout. (We will discuss this assumption later.) We 
will also consider the scalar potential along the real S and directions. Thus with 
s = \/2Re{S) and 0_ = v^2-Re($_), the tree-level scalar potential becomes 

„ M _) = ^ + ^-f) 2 . (3) 

The U(1)fi symmetry breaking transition occurs once s < s c = a/2^ 1//2 /A. 

In addition, we require the 1-loop effective potential as a function of s and 0_, 
AV(s, (/)-). This has been done previously for the case where 0_ is fixed at the 
minimum of its potential as a function of s However, as we wish to discuss the 
equations of motion of the s and <fi- field seperately here, we require the mass eigenval- 
ues and effective potential as a function of both fields. For s < s c the scalar and gauge 
boson mass terms and degrees of freedom as a function of s and 0_ are given in Table 1. 

lr The tree- level scalar potential for the D-term inflation model becomes equivalent to the minimal 
F-term inflation model 1 1 1| in the limit where A = \/2g. Therefore study of the D-term model should 
also give an insight into the growth of spatial perturbations in the F-term inflation model. 
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Table 1: Mass terms and physical degrees of freedom for s < s c . 



Here B = \\ 2 (s 2 + 2 ) + §<? 2 2 - g 2 i and C = |A 2 2 (§A 2 s 2 + §<? 2 2 _ - g 2 i 
|A 4 s 2 2 . r)± are the eigenstates of the si, <p\ mass matrix. (We define S = (si+is 2 )/V^ 
and = (0i + i<j> 2 ) / y/2.) A is the C/(1) F / gauge boson, with 2 the corresponding 
Goldstone boson 2 . 

The fermion mass eigenstates come from the mass matrix for \& = (x, A+, A_ , \ s ) 
where x is the U(1)fi gaugino, A + (A_) is the fermionic component of $+($_) and \ s 
is the fermionic component of the S chiral superfield, 
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The resulting eigenstates correspond to two Dirac fermons, of mass 
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where 
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X 2 



n 2 



2„2 j. 4 



(6) 



For s > s c there are two massive complex scalars <S>± and a single Dirac fermion 
\&, with squared masses given in Table 2. 



2 More precisely, $_ = -^e 10 , with the Goldstone boson degree of freedom. The mass terms of 
the physical scalar field 4> with 8 — are the same as those of <p\ with (f> 2 = 0. 
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Table 2: Mass terms and physical degrees of freedom for s > s c . 



The 1-loop effective potential as a function of s and is given by the Coleman- 
Weinberg formula 

^,0) = ^E(-i) F -^(^) , (7) 

where A is a renormalisation scale and the contribution is negative for fermions. An 
important question is how this should be interpreted when a scalar mass squared 
becomes negative for small field expectation values. This problem also arises in the 
case of the Standard Model Higgs effective potential. Weinberg and Wu JH] argue 
that only the real part of the effective potential corresponds to a correction to the 
energy density and so enters the equations of motion. Thus if m 2 < 0, we can simply 
substitute m 2 = |m 2 |e l7r in Eq. (jjj), in which case the real part of AV corresponds to 
replacing m 2 by |m 2 | throughout. We will use this proceedure in the following. 
The full equations of motion to 1-loop are then 

V 2 A 2 2 s dAV 

s + 3Hs -^ s = —2 aT> (8) 



and 

• V 2 ± AV0_ 2 / cj> 2 \ dAV 
0_ + 3H^ - — 0_ = ^- + g 2 f £ - ^ ] 0_ - ^- . (9) 



3 Spatial Perturbations and Initial Conditions 

3.1 Equivalent Classical Stochastic Field 

We are interested in the growth of spatial perturbations of the scalar fields. These have 
their origin in quantum fluctuations. In studying the growth of quantum fluctuations of 



a scalar field, the method used is that of the equivalent classical stochastic field (ECSF) 
[TK] . This is based on the observation that for quantum states which evolve to semi- 
classical final states, expectation values in the full quantum field theory are the same 
as those obtained by solving the classical equations of motion with initial conditions 
which have a classical stochastic distribution determined by the initial quantum state 3 . 
For the case we are considering here, the initial state is taken to correspond to the 
conformal vacuum state of a massless scalar field in (quasi) de Sitter space (H « 
constant). The corresponding initial conditions for the field modes are ^Hj 

i^ k )' 2= (iA?v)^(kyF = 2jb?( 1 + fcv) ■ 

Here we have used the notation of ^3], where y(k) is the Fourier transform of the 
perturbation of the scalar field </»(x), d>0(x) = y(x)/a (time dependence not explicitly 
stated), and we are considering initial conditions taken at scale factor a = 1. p(k) is 
the Fourier transform of p(x), the conjugate variable of y(x) 

p(x) = y-(fyy> ( 12 ) 

where a prime denotes differentiation with respect to conformal time rj (77 = / dt/a = 
— (aif) -1 ). These modes refer to Fourier transforms defined by 

y(x) = / y(k)e-^d 3 k , (13) 

T/-1/2 , 

vM = 7^m v vW e * x > (14) 

where we retain an explicit spatial volume V in order to account for the mass dimen- 
sions. For a real y(x) we have y(k) = y*{— k). 

Given the initial conditions, the subsequent evolution of the system is obtained 

by evolving the scalar field classically. In practice the important quantities are the 
3 A simple and useful illustration of this concept in the context of quantum mechanics is given in 

Hg. 



root mean square (r.m.s.) spatial averages of the field fluctuations. For example, once 
the r.m.s. fluctuation of the field amplitude is of the order of the homogeneous field 
amplitude, the field will almost entirely exist in the overdense regions and therefore we 
expect the spatial perturbations to go non-linear and to fragment to non-topological 
solitons. The mean squared fluctuation is given by 

(y(xf)=47rje\y(k)\ 2 f, (15) 

where k = |k|. Thus modes with |k| in a range Ak ~ k around k will contribute to 
the r.m.s. fluctuation squared an amount [1, 

A 2 ^(x) 2 )^4^ 3 |y(k)| 2 . (16) 

We will model the classical evolution of the r.m.s. fluctuation due to modes with 
|k| ~ k by considering a single mode with |k| = k and initial amplitude given by A&. 
We will be particularly interested in the dominant mode (the largest at late times), 
since its wavelength is likely to determine the size of the non-topological solitons 
which form once the inflaton sector energy density perturbations become non-linear 
[UJ |H] . (This is similar to what is observed numerically in the case of a complex scalar 
condensate and Q-ball formation ffiH I2T].) The modes of interest in the following will 
have wavelength small compared with the horizon at the end of inflation. In this case 
\kq\ = k/H > 1 and |y(k)| « 1/V2k. Then A fc « V^k. 

The momentum conjugate to y(x) in the limit k 3> H is given by 

p(x) = a(y(x)-ify(x)) , (17) 

where 'dot' denotes differentiation with respect to time. Thus if y(x)/y(x) = 50(x)/<50(x) > H 
(as expected if perturbation growth is significant) we have initially (at a = 1) p(x) ~ 
y(x) = <50(x). Therefore the initial mean squared rate of change of £0(x) correspond- 
ing to the range AA; ~ k around k is related to p(k) by (cf. Eq. ()16j0 . 

(50 2 (x)) fc « 47rfc 3 b(k)| 2 = 2nk 4 (|Y . (18) 
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Thus to model the ECSF initial conditions we will consider the evolution of a scalar 
field mode £0fc(x) with initially (at t = t ) 

6(j> k (x,t ) = A(t ) sin(k.x) , (19) 

where 

A(t ) = ^Phxk (20) 

and 

\A(t )\ = V^k 2 (I) . (21) 

(In practice we find that |A(i )| has a negligible effect.) 

We next consider the time at which the ECSF initial conditions should be applied. 
They strictly apply for the quantum vacuum state of an effectively massless scalar field 
(m <C H) in de Sitter space. For the inflaton, s, we may apply the initial condition 
at any time during slow rolling. So long as there is little growth of the effectively 
classical s fluctuation from this time until s reaches s c (which may be expected as 
there is only the 1-loop effective potential at this time) the same initial perturbation 
may be applied at s c (or more generally up to the time at which the U(1)fi symmetry 
breaking field develops a classical expectation value and the s perturbation begins 
to grow significantly). At s c the 0_ field is also massless. (However, there is a 4 
potential term which will limit the amplitude of effectively massless fluctuations to 
~ H 2 .) Thus a natural time to apply the equivalent classical initial conditions 
to both s and 0_ would be at s(t) ~ s c . 

At s c the mean value of the 0_ field is zero. The U(1) FI symmetry breaking occurs 
as a result of the growth of the 0_ quantum equivalent classical fluctuations when 
the 0_ field acquires a negative mass squared term once s < s c . The amplitude of 
the quantum equivalent classical fluctuation is <50_ q ~ \f2nk. We first note that 
the largest value of k for 0_ fluctuations for which the 0_ field can be considered 
as effectively massless corresponds to V" (5(f)- q ) (= 3g 2 5(f> 2 _ q /2) w H 2 . Thus the 
spectrum of equivalent classical fluctuations has an upper limit, k 2 < H 2 /37ig 2 , beyond 
which the quantum fluctuations will be suppressed by the effective mass term. A 
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fluctuation can grow due to the negative mass squared term so long as its physical 
momentum satisfies k 2 < \m^_\, where m^_ = (^f - — g 2 C)- I n addition, the mode 
can grow only once |m| _| > H 2 . Therefore, the largest amplitude mode which first 
begins to grow classically corresponds to k 2 ~ H 2 , which, for natural values of the 
U{1)fi gauge coupling not much smaller than 1, is of the same order as the upper 
limit of the spectrum of equivalent classical fluctuations. Thus we may consider the 
dominant mode to correspond to k ps H and to start growing once |m|_(s)| ~ H 2 , 
defined to occur at s — s±. Once the 0i interaction term is included, the minimum 
of the potential in the 0_ direction for a given s is at 4> m in = V / 2|^_|/fi f - Therefore 
the minimum of the potential is at around the amplitude of the largest quantum 
equivalent classical fluctuation which is able to grow due to the negative mass squared 
term. Once a mode is given a classical value of the order of (p m in, numerical solution 
of the classical equations of motion shows that a growing mode will remain close to 
4>min(s) as s decreases (Fig 1). This indicates that at any given time once s < s±, the 
dominant quantum equivalent classical 0_ mode will have amplitude around (p m i n (s). 

This also indicates that the size of the domains of U(Y)fi symmetry breaking and 
so seperation of the cosmic strings, which corresponds to the dominant 0_ fluctuation 
with k ~ H, will not be much smaller than the horizon (A ~ 2ii/k m if -1 ). In 
particular, the spacing of the cosmic strings will be much larger than the dynamical 
length scale of the inflaton sector fields, of the order of £ -1 / 2 . Therefore we expect 
that cosmic strings will be relatively rare and that they will not play a significant role 
in the subsequent evolution of the inflaton sector fields and energy density. 

3.2 Initial Conditions 

In order to analyse the growth of s perturbations in the presence of cosmic strings, it 
is necessary to perform a full lattice simulation. Here we wish to consider a simplified 
model without cosmic strings, which we have argued will not play a significant role 
in the growth of inflaton sector perturbations. This is achieved by correlating the 
direction of the initial 0_ quantum fluctuations. Thus we will consider an initial 
perturbation of the 0_ field with the same magnitude as the quantum equivalent 
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(a) g = 1.0 (b) g = 0.5 

Figure 1: After its initial perturbation (f> tracks its minimum (f> m i n {s) as a function of s 
as s decreases for both values of g. (Using j = 0.14 and k = 0.001s c .) 

classical fluctuations (i.e. of the order of m j„(s) at si) but homogeneous in space. We 
will refer to these as 'no cosmic string' (NCS) initial conditions in the following. The 
NCS initial conditions we consider are then 

s = si , (22) 
5s = ^Phxk ■ \5s\ = V2^k 2 (^j , (23) 

0- = <f>min(si) (24) 

and 

<50„ = O ; |<50_|=O, (25) 

with s 1 corresponding to \m^,(si)\ 2 rts H 2 . We have set <50_ = initially since the 
quantum equivalent classical 0_ perturbation has been rotated in space to become 
homogeneous. The initial value of the rate of rolling of the homogeneous inflaton field, 
s, which plays an important role in determining the subsequent growth of the spatial 
perturbations, is obtained by evolving the homogeneous s field from the slow-rolling 
regime at s ^> s c to s\ using the 1-loop effective potential. 
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The NCS initial conditions we have introduced are not strictly physical. Neverthe- 
less, the resulting growth of the inflaton perturbation has a general physical interpreta- 
tion, namely it is the maximum possible growth. This follows since the classical growth 
of inflaton perturbations is due to the non-zero 0_ expectation value as a function of 
s; in the limit 0_ = there is no significant growth since s then only experiences a 
1-loop effective potential as in the case s > s c . The effect of the cosmic strings is to 
reduce the value of 0_ at a given point in space relative to the case with no cosmic 
strings (the U{1)fi symmetry breaking field vanishing at the centre of the strings), so 
reducing the degree of U(l) fi symmetry breaking on average. Therefore we expect the 
growth of inflaton perturbations to be suppressed relative to the case with no cosmic 
strings. For the particular case of D-term inflation we have also argued that cosmic 
strings should be relatively rare on sub-horizon scale and so should not play a major 
role in determining the evolution of the energy density in the bulk between the strings. 
In this case the NCS initial conditions should be a good approximation to the fully 
realisitic case. 

4 Semi- Analytical Model for Spatial Perturbation 
Growth 

In this section we consider the growth of spatial perturbations using an analytical 
ansatz for the spatial dependence of the perturbation modes. Our method will be 
to study the evolution of a single perturbation mode of the s and 0_ field, using an 
ansatz for the functional form of the space dependence of the fields. We first discuss 
the ansatz for a generic real scalar field $. The field we consider is assumed to initially 
have the form $(x, t) = $(£) + 8Q(t) sin(k.x), where <£>(£) is the homogeneous field. 
In general, we would expect the spatial position of the extrema of $(x, t) at later 
times to remain the same as their initial positions. Thus we expect that $(x, t) = 
$(t)+/(t)g(x), where g(x) has the same periodicity in space as the initial perturbation. 
If the growth of the perturbation was strictly linear, then the correct choice would 
be g(x) = sin(k.x). However, it is possible that non-linearity could play a role in 
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the growth of the perturbation, in which case this assumption must be considered 
an ansatz, whose validity depends on how well the true interpolating function g(x) 
between the maxima and minima of $(x, t) in space can be modelled by sin(k.x). The 
match need not be exact for the physics of perturbation growth to be reasonably well 
modelled. 

The equation of motion of <3> is 



$ + -$ 

a 2 



dV($) 
d<S> 



(26) 



We introduce the ansatz $(x, t) = $(£) + 5&(t) sin(k.x), in which case the equations 
of motion become 

dV($) 



$ + 3H<f> + —5$ sin(k.x) 



(27) 



Let $ + and <3>_ be values of $(x, t) at points in space (sin(k.x) = ±1) corresponding 
initially to the maximum and minimum of $(x, t). Thus $ + = $(£) + 8Q(t) and 
$_ = $(£) — 5$(t). The equations at these points in space are then 



k 2 dV($) 



and 



$_ + 3#$_ - -^5$ 



<9$ 



ay ($) 



<9$ 



(28) 



(29) 



Thus with = ($ + — $_)/2 the equations for the evolution of the field at points in 
space corresponding to maxima and minima of $ are 



$ + + 3#$+ + —($+-$_) 



dV(&) 



and 



2a 2 



<&_) = 



9$ 



*4 



9$ 



(30) 



(31) 



Thus we have only to solve two coupled equations which are purely functions of t. We 
will refer to this as the semi-analytical model for perturbation growth. 

We can generalise the semi-analytical model to the case of hybrid inflation with 
two scalar fields. We consider an ansatz of the form s(x, t) = s(t) + 6s(t) sin(k.x) 
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and 0(x, t) = <fi(t) + 5<f)(t) sin(k.x). At s(t) = si the field at a point in space 
corresponding to the minimum of s(x, t), s_, will experience the largest symmetry 
breaking negative mass squared term. Therefore at this point 0(x, t) begins to grow 
first, corresponding to <f> + . (In accordance with NCS boundary conditions, we assume 
the symmetry breaking direction is correlated.) We therefore denote the fields at this 
point in space by (s_,0 + ), whilst the fields at the maximum of s(x, t) are (s + ,0_). 
(The spatial dependence of the s perturbation is therefore imprinted on the field.) 
This is illustrated in Figure 2. The semi-analytic equations are then 



s+ + 3Hs + + ^(s+ - s_) 



dV{sA) 



ds 



(32) 



and 



0- + 3#0_-— 



3 ^-"2^ (s+ - S - } 



k 2 

+ + — (0 + -0_) = 



dV{s, 



dv{sA) 



ds 
dV(s, 



(s+,<t>-) 



{s-,4>+) 



(s-,4>+) 



(33) 



(34) 



(35) 



Then 5s = (s + — s_)/2 and 50 = (0 + — 0-)/2, whilst the homogeneous mean fields 
are given by s(t) = (s + + s_)/2 and 0(t) = (0 + + <j)-)/2. 



5 Results for Perturbation Growth 

We next investigate the growth in inflaton perturbations using the NCS initial con- 
ditions and the semi-analytical ansatz for the evolution of the spatial perturbations. 
We will focus on the case g = 1 but include g = 0.5 for comparison. 

In the figures we show the evolution of the homogeneous mean field s(t) and the 
amplitude of spatial perturbation 5s(t) = (s + — s_)/2 for the semi-analytical model 
using the D-term hybrid inflation scalar potential with 1-loop corrections. (For clarity 
we show the magnitude of 5s(t) in the figures, such that 5s(t) is always positive.) 
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Figure 2: Schematic illustration of the time evolution of <j> and s as the inflaton s passes from 
slow-roll phase to U(1)fi symmetry breaking phase. 

The homogeneous inflaton field is evolved from the slow-rolling regime with initially 
s ps 20s c whilst the NCS boundary conditions for the perturbations are introduced 
at s\. We expect that the growth of the spatial perturbations will become fully non- 
linear with the formation of non-topological solitons ('fragmentation') once the spatial 
perturbation amplitude is similar to the homogeneous mean field, Ss(t)/s(t) ~ 1. 

In Figure 3 and 4 we show the evolution of Ss/s(t) for a mode of fixed wavenumber 
k = 0.01s c as X/g is increased for the case g — 1. (We consider a = 1 at s = s c , such 
that k is the physical wavenumber at this time.) In Figure 3a) with X/g = 0.1 there is 
little growth of 5s(t) whilst the mean field amplitude s(t) decreases due to expansion. 
In Figure 3b) with X/g = 0.12 significant growth of the spatial perturbation begins to 
be seen. In Figure 4a) with X/g = 0.14 we see that the spatial perturbation amplitude 
increases rapidly to become equal to the mean field, Ss(t)/s(t) ~ 1, in approximately 
one oscillation cycle of the homogeneous field. Finally, in Figure 4b) with X/g = 0.2 
the spatial perturbation grows more rapidly and dominates the homogeneous field at 
all times after about one oscillation cycle of the homogeneous field. 

These results indicate that the the growth of quantum fluctuations of the inflaton 
sector fields to form non-topological solitons can be rapid in realistic D-term inflation 
models, occuring in not much more than the time for a single coherent oscillation of the 
inflaton field. This is consistent with a tachyonic preheating picture of perturbation 
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growth [Z| rather than an inflaton condensate fragmentation interpretation 

We next consider the case of fixed X/g = 0.14 and investigate the effect of the 
variation of the mode k on the growth of the perturbations. Figure 5 shows the 
transition from a) the case of almost no-growth at k — 0.001s c to b) fully non-linear 
growth (5s/s(t) > 1.0) at k « 0.005s c . In Figure 6 from a) to b) a dramatic damping 
in the growth of the perturbation is seen as k increases above 0.01s c . Since the effect 
of large gradient terms in the equations of motion is to damp the growth of the 
perturbation, this would indicate that we have passed the maximum value, k ma x, for 
which growth of inflaton perturbations is possible. The dominant mode for a given 
X/g will correspond to that which grows to become fully non-linear first, which will 
correspond to k ~ k max . 

In Figures 7 and 8 we consider g = 0.5 whilst keeping X/g = 0.14. Again we vary k 
to show the transition from no growth to fully non-linear growth. In Figure 7b) we find 
that for k = 0.001s c the mode can grow to become non-linear but in a time much longer 
than the period of coherent oscillations. Growth of this mode would be well described 
by the inflaton condensate fragmentation approach jB] . However, it is the mode which 
grows most rapidly that will reach non-linearity first and so determine the length scale 
of the condensate lumps. Typically this dominant mode reaches non-linearity in a few 
oscillations. 

Comparing Figure 8a) with Figure 6a), we see that exactly the same maximum 
growth of the perturbation is obtained at k close to k max for both g = 0.5 and g = 1.0 
when X/g is fixed. This shows that the maximum growth of the perturbations is 
dependent only upon the ratio X/g. We also see that k max f& 0.007s c for g = 0.5, 
compared with k max ~ 0.015s c for g — 1, indicating that k max cx g for fixed X/g. 

An important quantity is the lower limit on X/g for which the perturbations can 
grow sufficiently to achieve fragmentation. We find that for X/g < 0.09 fragmentation 
does not occur for any value of k. Thus for X/g ^ 0.09 we expect fragmentation into 
non-topological solitons to occur in D-term inflation. The value of k max corresponding 
to the limiting case X/g ~ 0.09 is k max rs 0.008s c . 

Although the above results strictly apply to the case of D-term inflation, we note 
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Figure 3: k = 0.01s c . As we increase A from 0.10 to 0.12 we start to see growth in the 
perurbation amplitude though it is still well within the linear regime. 
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Fi gure 4: k = 0.01s c . (a) shows the critical value of — for which the perturbation just 
becomes fully non-linear. Any further increase ((b)) pushes the perturbation deeper into 
the non-linear regime. 
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that in the limit where the classical potential becomes identical to F-term inflation, 
A = V2g, condensate fragmentation is certain to occur in D-term inflation. This 
strongly suggests that fragmentation will be an inevitable feature of F-term inflation 
models. 



It is interesting to compare these results with the analytical results obtained in the 
inflaton condensate fragmentation approach 6j. Although the approximations made in 
|Hj are not valid for perturbations which grow rapidly on the time scale of the coherent 
oscillations, we nevertheless find that the dependence of k max and the fragmentation 
condition on A and g are well described by the analytical expressions in j^j. k max is 
given by 

V^A 2 

—s c , (36) 

where we have used R = s c and a Q = a in Equation 29 of 6j, consistent with k max at the 
onset of coherent oscillations. Thus we find k max cx g for fixed X/g. For X/g = 0.14 
we find from Eq. that k max = 0.008s c (0.016s c ) for g = 0.5(1.0), in reasonable 
agreement with our numerical results. The condition for fragmentation to occur from 
the inflaton condensate fragmentation analysis is 

g ~ M \l0) 

where (3 = log(s c /8s(t )). This depends only on the ratio X/g, in agreement with 
the numerical analysis. (The actual value of the lower bound from Eq. (|37j) can only 
be taken as suggestive, as the various approximations made introduce 'factor of 2' 
uncertainties.) We also note that the rate of growth of perturbation from the inflaton 
condensate fragmentation analysis (given by a (Equation 27) in 0) is proportional to 
k for values up to k max , beyond which there is no perturbation growth, in qualitative 
agreement with the numerical results. 

The broad agreement of the numerical k max with the inflaton condensate frag- 
mentation expression suggests that the exponential terms governing the growth of 
perturbations have the same parameter dependence in the realistic case as in the in- 
flation condensate fragmentation analysis, although there is no reason to expect the 
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perturbation growth at k ~ k max to be quantitatively well-described by this analysis. 

Finally, we would like to comment on the possible consequences of including the 
$+ field in the analysis. During inflation m$ + > \f2g£}l 2 ^> H. Therefore we expect 
that the homogeneous field will satisfy $ + = at the end of inflation. In this case the 
homogeneous <£>+ field will not evolve since $ + = is an extremum of the potential for 
all values of S and $_ . A quantum fluctuation of $ + could grow if the mass squared 
of $+ at $ + = were to become negative during the post-inflation evolution of the 
S and $_ fields. This could briefly occur as \S\ approaches zero if the value of the 
oscillating $_ field was also sufficiently larger than its minimum as a function of |*S'| at 
this time. However, the initial quantum de Sitter fluctuation of $ + will be suppressed 
by its large mass compared with H during inflation. In addition, the brief growth as 
\S\ passes through zero will result in only a small growth of the quantum fluctuation of 
$ + . This is particularly true since the dominant S modes reach fragmentation within 
only a few oscillations of the homogeneous S field. Therefore we expect that the <3> + 
field will have a negligible effect on the evolution of the energy density during the 
post-inflation era and that we may consider $ + = throughout. 
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Figure 5: ^ = 0.14, g = 1. As expected, for a very small value of k ((a)) there is no 
growth in the perturbation as the gradient term in the equations of motion are negligable. 
As k is increased ((b)) we get the expected evolution into the non-linear regime. 
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Figure 6: -| = 0.14, g = 1. Further increase ((a)) gives us the maximum value of k 
for which full non-linearity occurs earliest at this value of This corresponds to the 
dominant mode leading to fragmentation. Finally, (b) illustrates the case k ~ k max with 
suppressed perturbation growth. 
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Figure 7: ^ = 0.14, g = 0.5. As before we see the same increase in the rate of pertur- 
bation growth as we increase the initial mode size k. (b) shows a mode with slow growth 
towards full non-linearity. 
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Figure 8: ^ = 0.14, g = 0.5. (a) shows the limiting case with k « k max where the 
perturbation just reaches full non-linearity. Further slight increase in the mode size to 
k > k max ((b)) suppresses growth in the perturbation as before. 
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6 Consequences of Inflatonic Non- Topological Soli- 
ton Formation 



In this section we discuss some possible cosmological consequences of non-topological 
soliton (NTS) formation. After fragmentation the inflaton sector fields are expected to 
form NTS which we have labelled inflaton condensate lumps jSj, with essentially all of 
the energy density of the Universe concentrated in the lumps (HI |H] . Physically these 
are interpreted as a Bose condensate of scalars held together by an attractive self- 
interaction which balances against the gradient terms (uncertainty principle pressure) 
in the scalar field equations. The most obvious physical consequences of inflaton 
condensate lump formation will be related to the inhomogeneous distribution of the 
energy density of the Universe in the post-inflation NTS era and the inhomogeneous 
nature of reheating via NTS decay. 

The decay rate of the condensate lumps is related to (i) its stability as a solution of 
the classical equations of motion and (ii) the decay and annihilation rate of the scalars 
inside the NTS. For the case of a real scalar field it is likely that the condensate lumps 
are not stable, but instead slowly radiate scalar field waves until (depending on the 
potential) the amplitude of the scalar field inside the lumps becomes too small for 
the attractive interaction term to hold the lumps together. This has been studied 
extensively in [22], where the condensate lumps have been given the name 'oscillons'. 
The lifetime of the oscillons in a potential of the form V(4>) = m 2 2 /2 — A0 4 /4 (which 
approximates D-term inflation [U]) has been estimated to be around 10 3_4 m _1 (221 IE] 
In this case we would expect the scalars forming the lump to eventually disperse. 
However, during the period when the condensate lumps exist, the fact that the num- 
ber density of scalars inside the lumps is constant (unlike the case of a conventional 
homogeneous inflaton condensate jl]) means that perturbative inflaton annihilation 
and/or parametric resonance will be effectively enhanced and could dominate over 
perturbative inflaton decays as the main mode of reheating [Hj. This all depends on 
whether the classical decay of the condensate lumps is slow enough relative to annihila- 
tion/parametric resonance for the condensate lumps to have non-trivial consequences 
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for reheating. 

However, an important point is that in SUSY hybrid inflation models the inflaton 
is in fact a complex field. In D-term inflation there is a global U(l)s symmetry such 
that S — > e ia S and $ + — > e~ ta § + in the superpotential. As a result, there will be 
a Q-ball solution of the inflaton sector equations of motion [23 carrying a non-zero 
U(l)s charge. In this case it is likely that the initially neutral condensate lump will 
fragment into a ±Q-ball pair [21]. This occurs because although initially the inflaton 
field and condensate lump is real, small fluctuations of the phase of the inflaton field 
inside the lump can grow through mode-mode couplings, resulting in a seperation of 
charge in the quasi-stable condensate lump into a Q-ball, anti-Q-ball pair. This has 
been demonstrated in lattice simulations of a running mass chaotic inflation model 
[21] • Although the potential of this single-field inflation model is different from that of 
D-term inflation, we would expect the tendency of the system to reach its most stable 
configuration (±Q-ball pairs) in the presence of spatial perturbations to be repeated. 

Since the Q-balls are classically stable, reheating will be determined by the de- 
cay/annihilations of the globally charged inflaton sector fields forming the Q-balls i.e. 
reheating will occur via Q-ball decay, with the energy density of the Universe con- 
centrated inside the classically stable Q-balls until that time. Therefore reheating 
and thermalisation will be highly inhomogeneous in this case, unlike the case of a 
conventional homogeneous inflaton condensate. 

The formation of condensate lumps and fragmentation to Q-balls also results in 
a highly inhomogeneous energy density prior to reheating. One consequence of this 
is that the dynamics of scalar fields in SUSY and SUGRA cosmology will be altered. 
In SUGRA models, Planck-scale suppressed corrections result in a correction to the 
mass squared term of scalar fields the form \F\ 2 /M Pl « p/M Pl , where F represents the 
F-term of any field in the cosmological background and p is its the energy density. In 
the conventional case of a homogeneous coherently oscillating inflaton field following 
inflation, p is also homogeneous and so the correction to the mass squared term is of 
the form p/M Pl ~ cH 2 , where p the average energy density and |c| is of the order 
of 1 [Tlj. However, when the energy density is packed into Q-balls, the correction to 
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the mass squared term of the scalars outside the Q-balls will be effectively zero, whilst 
inside it will be much larger than H 2 . As a result, the dynamics of flat direction scalars 
in the MSSM [23] (and other so-called moduli fields) will be completely different from 
the conventional case. In the bulk outside the Q-balls (corresponding to most of the 
volume of the Universe) the scalar fields will not evolve until H 2 « m 2 susy (which occurs 
prior to reheating if the reheating temperature satisfies the thermal gravitino constraint 
Tr "Z 10 8 GeV where m 2 usy w (100 GeV) 2 is the conventional gravity-mediated 

SUSY-breaking mass squared term. Inside or close to the Q-balls, if we consider the 
mass squared correction due to p to be positive, the flat direction scalar will be driven 
towards zero. The full dynamics of the MSSM flat directions will require a complete 
solution of the scalar field equations in the presence of Q-balls, but we expect the 
dynamics in the bulk to be essentially that of c = in the conventional discussion 
of MSSM flat directions and moduli. This will have significant consequences for the 
cosmology of the MSSM, in particular for Affleck-Dine baryogenesis [231213 • We expect 
other consequences of the NTS dominated post-inflation era to become apparent in 
the future. 

7 Conclusions 

We have analysed the growth of spatial perturbations of the inflaton sector fields of 
D-term hybrid inflation including 1-loop radiative corrections, using a semi-analytic 
model with simplified initial conditions. We have argued that this model will give a 
good estimate of the growth of spatial perturbation modes in the fully realistic case. 

We find that the growth of the perturbations to fully non-linearity (indicating 
the onset of inflaton condensate lump formation) occurs once X/g ^ 0.09, and that 
the maximal growth factor of the perturbations is a function only of X/g. Therefore 
non-topological soliton formation is a natural possibility in D-term inflation models. 
The corresponding perturbation mode, which sets the length scale for the condensate 
lumps in the case where X/g ~ 0.09, is given by k « 0.008s c for g — 1, with k oc g for 
fixed X/g. The growth of spatial perturbations to full non-linearity is rapid, typically 
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after about one coherent oscillation of the inflaton field, which is consistent with a 
tachyonic preheating picture [7] of the growth of spatial perturbations. However, 
the value of k max , the largest wavenumber for which perturbation growth can occur, 
and the dependence of the fragmentation condition only on the ratio X/g are both 
consistent with the analytical predictions of the inflaton condensate fragmentation 
analysis 0. 

Although our analysis strictly applies to the case of D-term inflation, it strongly 
suggests that inflaton condensate fragmentation will be an unavoidable feature of F- 
term hybrid inflation, which has the same classical scalar potential as D-term inflation 
in the limit A = V2g. However, a full analysis including the 1-loop potential of F-term 
inflation remains to be done. 

The formation of non-topological solitons (inflaton condensate lumps) is undoubt- 
edly a feature of the cosmology of SUSY hybrid inflation models for a natural range 
of the inflaton sector couplings, resulting in a highly inhomogeneous post-inflation era 
with all the energy of the Universe packed inside the non-topological solitons. We have 
discussed some possible consequences of inflaton condensate lump formation for SUSY 
cosmology, including the formation of inflatonic Q-balls from decay of the condensate 
lumps, inhomogeneous reheating and the modification of flat direction scalar field dy- 
namics. We hope to develop a more detailed understanding of this post-inflation era 
in the future. 
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